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i \ In heavy-ion collisions, J ftp and Drell-Yan production are expected to be 

Ph. 
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c~| , affected by nuclear modifications to the free nucleon structure functions. If 

these modifications, known as shadowing, are proportional to the local nuclear 
density, the per nucleon cross sections will depend on centrality. Differences 
in quark and gluon shadowing will lead to a new source of impact parameter 
dependence of the J /ip to Drell-Yan production ratio. We calculate this ratio 
in the CERN NA50 acceptance with several shadowing parameterizations to 
explore its centrality dependence. 
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A significant 'anomalous' suppression of J /if) production has been observed in Pb+Pb 
collisions. The ratio of J /if) to Drell-Yan production is lower in central Pb+Pb collisions 
than extrapolations from more peripheral collisions, lighter ion interactions, and pA collisions 
suggest Centrality is inferred from the transverse energy, Et- 

Almost all calculations of J /if) and Drell-Yan production in nuclear collisions to date have 
been based on position independent nucleon |§ or nuclear structure functions ||. However, 
nuclear shadowing should depend on the parton's location inside the nucleus. If shadowing 
is due to gluon recombination [|J, nuclear binding or rescaling [g], or other local phenomena, 
it should be proportional to the local nuclear density. The only studies of spatial dependence 
have relied on qualitative measurements of impact parameter, such as dark tracks in emulsion 
||, to find evidence for such a spatial dependence. 

This letter presents calculations of the impact parameter dependence of nuclear shad- 
owing on J /if) and Drell-Yan production in heavy ion collisions. This spatial dependence 
has an important effect on the Et dependence of J / if> and Drell-Yan production. Since J / if) 
suppression is a predicted signature of quark-gluon plasma formation |7|, interpretations 
of the anomalous suppression in the NA50 Pb+Pb data |T[ should include the spatial 
dependence of the structure functions on the J /if} and Drell-Yan rates. We focus on the 
effects of shadowing and neglect J /if) absorption mechanisms. 

At leading order (LO), J /if) production is dominated by gluon fusion while Drell-Yan 
production is due to qq annihilation. Thus if the quark and gluon distributions are affected 
differently by shadowing, and if shadowing has a spatial dependence, the J /if) to Drell-Yan 
ratio can vary with impact parameter. We present results in the acceptance of the CERN 
NA50 experiment Jl|. 

To be consistent with the NA50 analysis, we calculate J / if) and Drell-Yan production at 
leading order, LO. The LO cross section for nuclei A and B colliding at impact parameter b 
with center of mass energy ^snn and producing a particle V (J/if> or 7*) with mass m at 
scale Q is 
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where ajj is the partonic ij — > V cross section. The nuclear parton densities, Ff-, are the 
product of momentum fraction, x, and Q 2 independent nuclear densities, pa\ position and 
atomic mass, A, independent nucleon parton densities, and a shadowing function, S l : 

F A (x, Q\ f, z) = p A {s)S\A, x, Q\ f, z)f, N (x, Q 2 ) , (2) 

where s = \/r 2 + z 2 . In the absence of shadowing, S i (A, x, Q 2 , f,z) = 1. The nuclear density 
is given by a Woods-Saxon distribution, pa(s) = po(l + lu(s/ Ra) 2 )/(1 + e^ s ~ RA '' d ) where 
electron scattering data || are used to fix Ra, d, u and p . 

We use the color evaporation model of J/ip production so that 



q=u,d,s 

,2 



The LO partonic cc cross sections are defined in |10| and m = XiX2$nn- The fraction of cc 
pairs that become J/^'s is fixed at next-to-leading order (NLO) ||. The ratio of the NLO 
to LO cross section is given by K t ^. We use the GRV LO [12| nucleon parton distributions 
with m c = 1.3 GeV and Q = m c as well as the MRS A' densities [[Tj| with m c = 1.2 GeV 
and Q = 2m c where m c and Q are chosen to agree with data ||. 

The LO Drell-Yan cross section depends on isospin since a™ ^ ^ a™ ^ 

f?M 2 )f?M 2 )ff- 2 = K e , v ^f- (4) 
x E 4[{%£(^ 2 ) + ^/^ , 

q=u,d,s 

where Za and Na are the number of protons and neutrons in the nucleus. We assume that 

fl = fd, fd = fu etc. In Eqs. (|) and (§, x lj2 = Qe**/y/s^ and Q = m - The factor 
accounts for the rate difference between the calculations and the data. 

Three parameterizations of shadowing, based on nuclear DIS [|TJ[], are used. The first, 
Si(A,x), treats quarks, antiquarks, and gluons equally without Q 2 evolution |Tj|. The 



other two evolve with Q 2 and conserve baryon number and total momentum. The second, 
Sl(A, x, Q 2 ), modifies the valence quarks, sea quarks and gluons separately for 2 < Q < 10 
GeV [|16[]. The most recent, S^(A, x, Q 2 ), evolves each parton distribution separately for 



Q > 1.5 GeV |T7|]. The S3 initial gluon distribution shows important antishadowing in the 
region 0.1 < x < 0.3 with sea quark shadowing in the same x range. In contrast, S2 has less 
gluon antishadowing and essentially no sea quark effect. 

Since we assume that shadowing is proportional to the local nuclear density, the spatial 
dependence is defined as 

S* ws = S'iA, x, Q 2 , f, z) = l + N WS [S\A, x, Q 2 ) - 1]^ , (5) 

Po 



where iVws = is chosen so that (1/A) J d 3 sp(s)Sws — S\ similar to flnj. For lead, Nws 



1.32. At large radii, s 3> Ra, S\ys ~~ * 1 while at the nuclear center, the modifications 
are larger than the average S l . An alternative parameterization, S^, proportional to the 
thickness of a spherical nucleus ||19| , leads to a slightly larger modification in the nuclear 
core. 

We model our calculations to the NA50 experimental acceptance with center of mass 
rapidity < y cm < 1 and decay angle in the Collins-Soper frame | coshes I < 0.5 [[!]]. The 
Drell-Yan spectrum is measured for m > 4.2 GeV and the factor K exp is extracted by 
comparing to Eq. (H) calculated with the GRV LO distributions flTJ. The ratio a J ^/a DY 
is formed by extrapolating the calculations to 2.9 < m < 4.5 GeV with the same K cxp 
because J/ip and ip' decays dominate the region 2.7 < m < 3.5 GeV. Table | give the 
impact parameter averaged cross sections per nucleon pair for the GRV LO and MRS A' 
distributions, for pp, Pb+Pb and S+U interactions. The numbers show the effects of isospin 
and shadowing at y/s^jy = 17.3 GeV. With the MRS A' set, the isospin correction is quite 



small in the extrapolated region | 2T| . Because isospin is unimportant in J/ip production, 
only the Pb+Pb cross section is shown. However, the J ftp data suggests that the AB 
dependence might be stronger if shadowing could be removed from the data. Although 
the choice of parton densities influences the isospin effect and K exp , the average Drell-Yan 



shadowing is changed by less than 1% while the difference in J/ip shadowing can be as large 
as 5%, largely because the gluon distribution is imprecisely measured. The table shows that 
the dependence on the nuclear species is weak. The change in the shadowing is also ~ 1% 
when ^snn = 19.4 GeV. Most important is the mass interval: shadowing is 5% stronger in 
the measured region than the extrapolated region. 

To illustrate how shadowing could affect the interpretation of J/ip suppression, Fig. [I] 
shows the Drell-Yan mass distribution for m > 2.9 GeV in three impact parameter bins 
relative to no shadowing, S = 1. Shadowing changes the slope of the spectrum, producing 
a ~ 20% change in the predicted rate at m ~ 9 GeV. The absence of Q 2 evolution causes 
the S\ results to decrease faster with mass than S2 or 5*3. Including only spatially-averaged 
shadowing increases K exp over that needed for S = 1 in the measured region relative to 
the extrapolated region, shown in Table |, as well as further increase the discrepancy in 
central collisions while overestimating K exp in peripheral collisions. At small radii, Fig. |l](a), 
da BY I dm drops more rapidly than the impact parameter averaged spectra. When b ~ Ra, 
Fig. |l|(b), the averaged and spatial dependent spectra approximately coincide. At large 
radii, Fig. [TJ(c) , shadowing is reduced, approaching the unshadowed spectra. These results 
show that using a calculation to extrapolate to an unmeasured region is problematic. 

Figure |2| shows the impact parameter dependence of the J / ip and Drell-Yan cross sections. 
The X2 ranges of the two processes nearly coincide: the Drell-Yan region 2.9 < m < 4.5 GeV 
corresponds to 0.062 < X2 < 0.26 while the J/ip cover 0.066 < x?, < 0.18. Since all three 
parameterizations assume some gluon antishadowing in this region, a J ^ is always enhanced 
in Pb+Pb collisions. Because Si is the same for quarks and gluons, J/ip and Drell-Yan are 
equally affected by shadowing. On other hand, £2,3 < 1> reducing a DY . 

In Fig. U, a J ^/a DY , calculated in Eqs. is presented as a function of Et- The 

correlation between Et and b is based on the number of nucleon participants , in agreement 
with the most recent NA50 Et distributions [p0[] . The Drell-Yan cross section is corrected 
for isospin and a J ^ /a DY is calculated at ^snn = 19.4 GeV. After these adjustments, with 
5=1, cr J /^/cr DY ~ 40.3, in agreement with the published NA50 data The combined 
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Drell-Yan shadowing and J/ip antishadowing in Fig. [| increases cr J ^ /<r DY to 40.6 with Si, 
44.5 for S 2 and 54.4 using S3. The Si ratio is independent of E T because the b dependence 
cancels. However, S2 and 5*3 vary with Et- The S2 ratio rises about 7% while the 53 ratio 
increases ~ 11% as (Et) grows from 14 GeV to 120 GeV. These enhancements are opposite 
to the observed drop at large Et and small b neglecting shadowing. 

Because of uncertainties in the gluon shadowing parameterization, it is difficult to draw 
detailed conclusions. However, Si should represent a lower limit and 5*3 an upper limit. 
A stronger spatial dependence such as would slightly increase the effect with Et while 
parameterizations based on e.g. nuclear binding || might predict a smaller effect. 

In conclusion, we have studied the impact parameter dependence of the ratio a J ^ /cr DY 
using a spatially dependent shadowing model. We find that the ratio increases at small b 
(large Et) compared with more peripheral collisions. The magnitude of the effect depends 
on the chosen parameterization. Neglecting shadowing could lead to an increased K exp at 
y / J/v r /v = 17.3 GeV since the measured cross section is more strongly affected by shadowing 
than the extrapolated cross section. In addition, using an impact parameter averaged spec- 
tra in central collisions would tend to underestimate the total number of Drell-Yan pairs, 
increasing a J ^ /a DY . If this effect could be identified and corrected for in the data, then 
cr J// ^/cr DY would rise ~ 10% at low Et and drop ~ 4% at high Et, enhancing the discrep- 
ancy between absorption models and the data. At higher y/s^N, such as at future heavy-ion 
colliders, the shadowing effect will be larger [21] since these colliders probe lower x values. 
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TABLES 



Drell-Yan (2.9 < m < 4.5 GeV) Drell-Yan (4.2 < m < 9 GeV) J/if> 

o pp (pb) apbPb (pb) asu (pb) Opp (pb) cr PbPb (pb) cr S u (pb) cr Pb pb (nb) 











GRV LO 








S = l 


16.6 


12.7 


13.3 


2.31 


1.66 


1.76 


1.61 


s = s 1 


— 


12.8 


13.5 


— 


1.61 


1.72 


1.64 


s = s 2 


— 


12.5 


13.1 


— 


1.56 


1.67 


1.84 


s = s 3 




11.9 


12.6 




1.48 


1.61 


2.04 










MRS A' 








S = l 


18.8 


18.3 


19.0 


2.15 


2.17 


2.29 


1.53 


S = Si 




18.6 


19.4 




2.10 


2.26 


1.58 


s = s 2 




18.1 


18.9 




2.03 


2.19 


1.74 


s = s 3 




17.2 


18.1 




1.93 


2.10 


1.85 



TABLE I. The Drell-Yan extrapolated and measured cross sections in pp, Pb+Pb and S+U 
collisions and the J/ip cross section in Pb+Pb collisions in the NA50 acceptance with the GRV LO 



and MRS A' parton densities at ^/snn = 17.3 GeV. We have not included K exp = 2.4 for GRV 
LO and 1.7 for MRS A'. 
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FIG. 1. The Drell-Yan dilepton mass spectrum, da BY /dM, without and with spatial depen- 
dence relative to S = 1. The curves correspond to 6-averaged results with S\ (dashed), S2 (dotted) 
and 53 (dot-dashed). The spatial dependence is illustrated for Si 5 ws (circles) S2,ws (squares) and 
^3,WS (triangles). The impact parameter ranges are (a) < b < 0.2Ra fm, (b) 0.9Ra < b < 1.1 R a 
fm and (c) 1.9Ra < b < 2.1Ra fm respectively. 
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FIG. 2. The (a) J/ip and (b) Drell-Yan rates in Pb+Pb collisions at ^/sj^ = 17.3 GeV 
relative to production without shadowing, S = 1, as a function of impact parameter. The curves 
and symbols are defined in Fig. [I]. 
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FIG. 3. The ratio of J/ip to Drell-Yan production, as a function of transverse energy, Et- The 
curves and symbols are defined in Fig. [l|. 
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